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ПРО АНАЛІТИЧНУ АПРОКСИМАЦІЮ ГУСТИНИ 
ЗА ДАНИМИ ГРАДІЄНТА СИЛИ ТЯЖІННЯ 


Математичні постановки конкретних геофізичних задач зводяться у підсумку до визначення зна- 
чень коефіцієнтів рівнянь математичної фізики, які входять у ці постановки. В руслі цього факту розв'я- 
зання задачі Алексідзе у вигляді нелінійного інтегрального рівняння зведено до лінійної комбінації шу- 
каних розв'язків. Ця альтернативна задача зведена до класичної задачі варіаційного числення. 


Ключові слова: гравіметрія; обернена задача; аналітична апроксимація; аналітичне продовження; 
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Обмеження математичного моделювання в 
геофізиці. Фахівці функціонально-аналітичного 
напряму геофізики розробляють методи і алгорит- 
ми в рамках відомих інтегро-диференціальних рів- 
нянь математичної фізики, які описують поведінку 
геофізичних полів. Їх зусилля спрямовані на виз- 
начення деяких коефіцієнтів цих рівнянь та внесе- 
ння відповідних поправок за "складність" геологі- 
чного середовища. Ці коефіцієнти визначаються 
за допомогою ітераційних схем в рамках методики 
розв'язання некоректних обернених задач геофізи- 
ки за Тихоновим |Лаврентьев, 2010). 

Проте визначення цих коефіцієнтів за функція- 
ми, які утворюють відповідні системи рівнянь - 
принципово неоднозначне. Це очевидно при сері- 
йних машинних розрахунках на системах великої 
розмірності | Страхов, 2007|. А при найменших по- 
хибках у вхідних даних розв'язки ще Й істотно не- 
стійкі. Здолати цю нестійкість на порядок важче, 
ніж отримати формальні розв'язки задач. До того 
ж, слід вхідних функцій, за яким визначають їх кое- 
фіцієнти - двовимірний, а відновити потрібно три- 
вимірне поле. А вже за значеннями поля слід знай- 
ти коефіцієнти при третій координаті, яка зовсім 
не співпадає з розмірністю сліду. 

Формальний розв'язок цієї проблеми дає псев- 
допросторовий розподіл коефіцієнтів (густин чи 
інших параметрів), який "розчиняється" серед ек- 
вівалентних розв'язків задачі. Крім того, позичив- 
ши апарат математичної фізики, геофізики одно- 
часно успадкували з ним і евклідову структуру 
просторів функцій, яка для регіональних побудов 
недостатня |Дубовенко, 2011а|. Але обмеженість 
інтегро-диференціальних операторів, які апрокси- 
мують математичну модель геологічного середо- 
вища - лише методична сторона проблеми. 

Вхідні дані - це, загалом, двовимірні проекції 
геофізичних полів на поверхню Землі, виміряні з 
похибками на розрідженій нерегулярній мережі 
спостережень. Вони апріорі отримані з невідоми- 
ми похибками, тому на додачу до вказаних вище 
методичних проблем ускладнюють поведінку (1 
визначення) функцій, які входять у відповідні мо- 
дельні рівняння. А розвиток комп'ютерного моде- 
лювання призвів до застосування дискретної шка- 
ли вхідних функцій. Враховуючи, що досі не ос- 
воєно пряме отримання вхідних даних для оберне- 
них задач геофізики із баз даних пунктів спостере- 


жень |Якимчик, 2010|, переведення неперервних 
вхідних даних у дискретну форму в цифровій кар- 
тографії неухильно вносить додаткові спотворен- 
ня, крім похибок вимірів. Ці проблеми являються 
наслідком методології отримання даних минулого 
століття, яка відживає своє. Ми їх називаємо мето- 
дологічними проблемами. 

Новий  геофізико-математичний апарат. 
Указані вище, та інші математичні і методичні 
проблеми, зумовлені неадекватністю апарата ма- 
тематичного моделювання геофізичних даних, 
спонукають геофізиків розробляти власний геофі- 
зико-математичний апарат. В його рамках істотно 
перероблені відомі математичні методи з урахува- 
нням фактів, які властиві винятково геофізичним 
явищам. Розпочав розробку своєрідного "геофізи- 
чного діалекту" В.М. Страхов |Страхов, 2001) у 
сфері потенціальних полів, врахувавши ортогона- 
льність сигналу і похибки та дискретність даних; 
продовжив В.В. Гольдін в сейсміці, В.В. Аксьонов 
у електродинаміці і т.д. 

Без чіткої національної програми розвитку тех- 
нологій потенціальних полів і створення мережі 
моніторингу цих полів не можна впливати на ре- 
форму принципів отримання вихідних даних. Мо- 
жливо, адекватна потребам геофізиків-польовиків 
та інтерпретаторів система неперервних спостере- 
жень для формування загальнодоступної бази да- 
них гравімагнітних полів постане в ближчій перс- 
пективі. Але можна подати деякі міркування щодо 
нових зображень математичних моделей геологіч- 
ного середовища і геофізичних полів у рамках ап- 
роксимаційного підходу |Страхов, 2001) до ство- 
рення подібних конструкцій. 

Нова постановка оберненої задачі гравімет- 
рії. Стійкі способи відновлення потенціалу сили 
тяжіння вкрай важливі 1 для детальної гравіметрії і 
для глобального відновлення фігури Землі. Остан- 
нім часом ради підвищення стійкості та надійності 
розв'язків у практичних застосуваннях задач тео- 
рії потенціалу вживають нову величину - градієнт 
потенціалу сили тяжіння. Теоретично вона визна- 
чена ще в роботі | Алексидзе, 1965|. Сучасні трак- 
тування градієнта потенціалу сили тяжіння орієн- 
товані на застосування у поєднанні з класичним 
зображенням аномалії сили тяжіння як вертикаль- 
ної похідної потенціалу. 

Одначе значно вищу стійкість і ширше коло за- 


стосувань (хоча й складніше двоступеневе обчис- 
лення) має аналітичне продовження згаданого гра- 
дієнта для спеціально отриманої |Дубовенко, 20- 
11а| нової аналітичної (диференціальної) моделі 
поля сили тяжіння. В ній врахована векторна при- 
рода сили тяжіння, градієнт її зміни залежно від 
рівня кривизни поверхні вимірів. Також визначені 
похибки, які виникають при заміні прямого опера- 
тора відомого в аналітичному продовженні рівня- 
ння Пуассона на рівняння сили тяжіння: вони про- 
порційні геометричним параметрам області дослі- 
джень. 

Коротко основні етапи її виведення ілюструє 
такий ланцюжок рівнянь: 
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В рамках розв'язання однорідної частини цьо- 
го рівняння виникла нелінійна гранична задача 
Алексідзе для рівняння Лапласа |Дубовенко, 20- 
116), яка містить в граничних умовах значення мо- 
дуля градієнта потенціалу сили тяжіння. Її розв'я- 
зок повністю теоретично обгрунтований для тяжі- 
ючих тіл, обмежених поверхнями Ляпунова. Певні 
труднощі її розв'язання відомими чисельними ме- 
тодами виявлені у роботі | Дубовенко, 2010). 

Загалом, для визначення модуля градієнта сили 
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вираз (2), визначається з нелінійного рівняння: 
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Методика розв'язання задачі. Щоб отримати 
розв'язок рівняння (3), в |Дубовенко, 2011а| зап- 
ропоновано дві ітераційні схеми, на жаль, мало 
придатні для практичних обчислень великої розмі- 
рності. Гадаємо, це сталось внаслідок слабкої збі- 
жності відповідного інтеграла (в смислі головного 
значення). А його регуляризація через зображення 
відрізком ряду не гарантує бажаної точності роз- 
в'язків. 

Один із можливих способів розв'язання рівня- 
ння (3) такий. Подамо його шуканий розв'язок 
с5), 5є ОС у вигляді такої лінійної комбінації 
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В такому поданні задача в підсумку зводиться 
до класичної апроксимаційної постановки: знайти 


коефіцієнти аї рівняння (4) так, щоб виконува- 
лась умова мінімізації нев'язки функціонала 
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Висновок. Тепер, зважаючи, яку з популярних 
у теорії обернених задач метричних норм ми обе- 
ремо, задача знаходження аналітичної апроксима- 
ції (4)-(5) вирішується у тому чи іншому просторі 
за допомогою класичних варіаційних методів, зок- 
рема, модифікацій метода Ньютона. Розв'язок от- 
римується у кілька разів швидше (залежно від ар- 
хітектури комп'ютера), а точні переваги виявлять- 
ся після випробувань на польових даних. 
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ОБ АНАЛИТИЧЕСКОЙ АППРОКСИМАЦИИ ПЛОТНОСТИ 
ПО ДАННЬІМ ГРАДИКЕНТА СИЛЬІЇ ТЯЖЕСТИ 


ЮХ.Й. Дубовенко 


Математические постановки конкретньвіх геофизических задач сводятся в итоге к определению 
значений козффициентов уравнений математической физики, входящих в зти постановки. В русле зтого 
факта решение задачи Алексидзе в виде нелинейного интегрального уравнения сведено к линейной ком- 
бинации ИСКОМЬХ решений. Зта альтернативная задача сведена к классической задаче вариационного ис- 
числения. 


Ключевье слова: гравиметрия; обратная задача; аналитическая аппроксимация; аналитическоє 
продолжениєе; градиент силь тяжести; линейная контактная задача; задача Алексидзе; метод Ньютона. 
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